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Abstract 

A fascinating inverse problem that has been receiving considerable attention is the 
construction of realizations of random two-phase heterogeneous media with a target 
two-point correlation function. However, not every hypothetical two-point correlation 
function corresponds to a realizable two-phase medium. Here we collect all of the known 
necessary conditions on the two-point correlation functions scattered throughout a di- 
verse literature and derive a new but simple positivity condition. We apply the necessary 
conditions to test the realizability of certain classes of proposed correlation functions. 

1 Introduction 

Random two-phase heterogeneous media abound in synthetic products and nature. Exam- 
ples include composite materials, colloidal dispersions, gels, foams, wood, geologic media, 
and animal and plant tissue.^ 12 El IH 13 The effective transport, mechanical and electro- 
magnetic properties of such heteroegeneous materials are known to depend on correlation 
functions that statistically characterize the microstructure.j^ It has recently been suggested 
that microstructure reconstruction problems can be posed as optimization problems. jHl E] A 
set of target correlation functions are prescribed based upon experiments or theoretical models. 
Starting from some initial realization of the random two-phase medium, the reconstruction 
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method proceeds to find a realization by evolving the microstructure such that the calculated 
correlation functions best match the target functions. This intriguing inverse problem is solved 
by minimizing an error based upon the distance between the target and calculated correlation 
functions. The two-phase medium can be a dispersion of particles in some matrix (liquid or 
solid) 1^ or, more generally, a digitized image of a two-phase material. [7| 

An effective reconstruction procedure enables one to generate accurate structures at will, 
and subsequent analysis can be performed on the image to obtain desired macroscopic prop- 
erties (e.g., transport, electromagnetic and mechanical properties) of the media. This be- 
comes especially useful in generating three-dimensional structures from planar information 
when three-dimensional imaging techniques are not available: a "poor man's" tomography 
experiment. [HI E] 

Interestingly, the same procedure has been used to "construct" realizations of two-phase 
media from a hypothetical target correlation function. [71131101121 In this mode, the procedure is 
referred to as a construction algorithm. There are mnay different types of statistical descriptors 
of two-phase media, [3^ but the most basic one is the two-point correlation function, which 
gives the probability of finding two points in one of the phases (see definition below) and is 
obtainable from small-angle X-ray scattering. ^I] The construction algorithm can be employed 
to determine if a prescribed two-point correlation function is in fact realizable. If such a two- 
point correlation function is realizable, then the procedure could be used to categorize classes 
of random microstructures, which would be a valuable accomplishment. However, not every 
hypothetical two-point correlation function corresponds to a realizable two-phase medium [S] . 
Therefore, it is of great fundamental and practical importance to determine the necessary 
conditions that realizable two-point correlation functions must possess. [121 IS] We note in 
passing that this question is closely related to realizability issues of pair-correlation functions 
of many-particle systems. jl31 IT^ [T3] 

One aim of this paper is to gather all of the known necessary conditions on the two-point 
correlation function of two-phase random media, also known as "random closed sets" in the 
field of stochastic geometry. [T7j Some of these conditions are well-known in the physical 
sciences literature, but others are more arcane and are contained in obscure mathematical 
technical reports and/or proceedings. We also derive a new but simple positivity condition 
on the two-point correlation function. We consider some illustrative examples of proposed 
correlation functions and test whether they can correspond to realizable two-phase random 
media. 
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2 Necessary Conditions 



Here we collect all of the known necessary conditions on the two-point correlation function 
of random media that are scattered throughout a diverse literature. We also derive a new 
positivity condition. 

Each realization uj of the two-phase random medium occupies the region of c?- dimensional 
Euclidean space V G S?'^ of volume V that is partitioned into two random sets or phases whose 
interiors are disjoint: phase 1, a region Vi{uj) of volume fraction 0i, and phase 2, a region 
V2{^jj) of volume fraction 02- For a given realization the indicator function X^^\x; u) for 
phase i at any position vector x G V is defined by 



Thus, a two-phase random medium is described by a binary stochastic process {T'-*)(x) : x G 
3fJ*^}. For statistically homogeneous but anisotropic media, the first two correlation functions 
are given by[S] 



where i = 1 or 2, angular brackets denote an ensemble average, r = xi — X2, and the symbol 
u is henceforth dropped for brevity. (The generalization to n-point correlation functions for 
n > 1 is straightforward. Clearly, (pi lies in the closed interval [0, 1] and 0i + 02 = 1- The 
two-point or autocorrelation function 5*2* ""(r) for statistically homogeneous media gives the 
probability of finding the end points of a vector r in phase i. Debye and Bueche showed 
that the two-point correlation function of a porous medium can be obtained experimentally 
via small X-ray scattering. Note that the two-point function for phase 2 is simply related to 
the corresponding function for phase 1 via the expression 
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and 




(3) 




(4) 



and thus the autocovariance function 




(5) 



for phase 1 is equal to that for phase 2. Generally, for r = 0, 




(6) 
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and in the absence of no long-range order, 



lim Si'\r)^<pl (7) 

|r|— >oo 

An important necessary condition for the existence of a two-point correlation function 
5*2* (r) for a two-phase statistically homogeneous medium d dimensions is that the (i- dimensional 
Fourier transform of autocovariance xl^")) denoted by must be nonnegative for all wave 

vectors, 12] i.e., 

X(k) = / xir)e~^^ dr>0, for all k, (8) 

where x(r) is given by (jSJ. Physically, this nonnegativity condition arises because is 
proportional to the scattered intensity, which must be positive. |3j This is sometimes called 
the Wiener-Khint chine condition, which is necessary but not sufficient for the class B of 
correlation functions that come from binary stochastic processes {X(*)(x) : x G 9?*^}. The 
Wiener-Khintchine condition is easily proved by exploiting a well-known theorem that states 
any continuous function x(r) must be positive semi-definite in the sense that for any finite 
number of spatial locations ri, r2, . . . , in 3?'^ and arbitrary real numbers Oi, 02, . . . , a^, 

m m 

Y,Y.^i^jX{ri-rj) >0 (9) 

i=i j=i 

if and only if it has a nonnegative Fourier transform Note that this property does 

not prevent from being pointwise negative for certain r. Importantly, whereas the real- 
space condition is difficult to check, the spectral version (jH)) is straightforward to test. It 
is noteworthy that if the medium in d dimensions is both statistically homogeneous and 
isotropic, then the one-, two-, ■ ■ ■ and (i-dimensional Fourier transforms of x(r) niust all be 
nonnegative. |12j This is a consequence of the fact that x(r) for such a random medium is an 
invariant in any m-dimensional subspace, where m = 1,2, . . . ,d — 1. 

The task of determining the necessary and sufficient conditions that B must possess is 
very complex. It has been shown that autocovariance functions in B must not only meet the 
condition of (jH)) but another condition on "corner-positive" matrices. [3] Since little is known 
about corner-positive matrices, this theorem is very difficult to apply in practice. Thus, a 
meaningful characterization of B remains an open and interesting problem, especially in the 
context of (i-dimensional two-phase random media. 

No attempt will be made to address the complete characterization of B here but instead 
we summarize the known necessary conditions, in addition to condition (jS)), that characterize 
B, most of which are described in Ref. (3). The two-point correlation function must satisfy 
the bounds 

< S^\r) < (f), for all r. (10) 
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The lower bound states that 5*2 (r) must be nonnegative for all r, but we show below that 
either S'2^''(r) or S^\r) must strictly be positive for 0j ^ 1/2. The corresponding bounds on 
the autocovariance function is given by 

— min(0^,02) ^ < 0i02 for all r. (11) 

Another consequence of the binary nature of the process in the case of statistically homoge- 
neous and isotropic media, i.e., when 5*2*^ (r) only depends on the distance r = |r|, is that its 
derivative at r = is strictly negative or 

< for all < 0i < 1. (12) 

r=0 

This is a consequence of the fact that slope at r = is proportional to the negative of the 
specific surface. |3] This means that 5*2* ""(r) has a cusp at the origin, implying that the two- 
point function is nonanalytic at the origin. It is a property of binary processes that if Xi(r) 
and X2(r) are in B, then Xi(r) ■ X2(i') ^ B and a;Xi(r) + (1 — «)X2(i") G B for every a G [0, 1]. 
This was proved by Shepp^HI in one dimension, but the proof should extend trivially to d 
dimensions. 

A little known necessary condition for statistically homogeneous media is the so-called 
"triangular inequality" first derived by Shepp [TH] and later rediscovered by Matheron: [12] 

^«(r)>S«(s) + ^«(t)-0„ (13) 

where r = t — s. The derivation of the triangular inequality (fT^ is straightforward. Following 
Shepp, we introduce the random variable 

/■^ f 1; if X G Vj, 

fW(x) = 2jW(x) - 1 = i ■ (14) 

[ — 1, otherwise. 

The mean of y(^)(x) is (F(^)(x)) = 20i - 1, which is equal to zero if 01 = 02 = 1/2- Observe 
that y«(xi) -F« (X2) + y*^*^(x3) is an odd number (either —3, —1, 1 or 3) and therefore 

([r«(xi) - F«(X2) + F«(X3)]2) > 1. (15) 

Using the fact that (y'^*)(xi)y^*)(x2)) = 4S'2*''(xi — X2) — 40j + 1, where we have invoked 
statistical homogeneity, we immediately obtain the triangular inequality (fTSj) . 

Note that if the autocovariance xi'^) of a statistically homogeneous and isotropic medium 
is monotonically decreasing, nonnegative and convex (i.e., d'^x{^)/dr'^ > 0), then it satisfies 
the triangular inequality (jl3p . The triangular inequality implies a number of pointwise 
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conditions on the two-point correlation function. For example, for statistically homogeneous 
and isotropic media, the triangular inequality implies condition (|T2|l . the fact that the steepest 
descent of the two-point correlation function occurs at the origin, [THj i.e., 

\S2\0)\ > \si'\r)\ for all r, (16) 

and the fact that 5*2*'' (r) must convex at the origin. pU] i.e.. 

From the "stochastic continuity" theorem for general stochastic processes, |2L it follows that 
if 5*2*^ (r) is continuous at r = 0, then it is continuous for all r. This continuity condition can 
also be proved using the triangular inequality. Note that 5*2*^ (r) can be discontinuous at the 
origin if the specific surface s is infinitely large. 

The triangular inequality is actually a special case of the more general condition [TS] 

m m 

fjil, i = l,...,m, m odd. (18) 

i=i j=i 

This necessary condition is much stronger than ©, implying that there are other necessary 
conditions beyond the ones identified so far. However, the condition ()18|) is difficult to check 
in practice because it does not have a simple spectral analog in contrast to (jHl) [cf. (jH))]- Note 
that the integers = ±1 in (jl8j) can be replaced with general integers, which would lead to 
an even more general condition on x(r). 

Here we report a new simple consequence of the lower bound of expression ()11|). Because 
the autocovariance is the same for phase 1 and phase 2, then it immediately follows from the 
lower bound of (fTT|) that 

S^\r) > max(0, 20^ - 1) for all r. (19) 

Thus, for 0j > 1/2, 5*2*^ (r) is strictly positive such that it must be greater than 20j — 1. 
Interestingly, the lower bound of (fTT|) for the autocovariance x(r), first obtained in Ref. 3, 
was derived from the trivial pointwise nonnegativity condition S'2*'*(r) > 0. However, the 
consequences of going back to the two-point correlation function S'2*''(r) were heretofore not 
examined. The nontrivial positivity condition (jl9p arises because the statistics of phase 1 
are not independent of the statistics of phase 2. Since 0^ is the large-distance asymptotic 
limit of S'2*''(r), its globally minimum value or, more precisely, its infimum (greatest lower 
bound) must be less than or equal to 0?. (Technically, one must consider the infimum and 



> 0. 



(17) 
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not the minimum because the minimum may not actually be achieved, e.g., a monotonically 
decreasing function that only asymptotically approaches its minimum value of (pf.) Clearly, 
the lower bound holds for the infimum of 5*2*^ (r), which will be denoted by inf [5'2*^(r)]. 
In summary, the infimum of any two-point correlation function of a statistically homogeneous 
medium must satisfy the inequalities 

max(0, 2(j)i - 1) < M[S^^\r)] < 0^ (20) 

(see Figure 1). 




Figure 1: Graphs of the upper bound (dashed curve) and lower bound (solid lines) of on 
the infimum of 5*2*'* (r) for a statistically homogeneous medium. 



3 Illustrative Examples 

It is convenient to introduce the scaled autocovariance function /(r) defined as 

/(r) = -— — = — for < r < +oo. 

0102 0102 

From p9|) . we obtain the triangular inequality for / to be 



:2ii 



/(r)>/(s) + /(t)-l. 



(22) 



Moreover, the bounds become 



— mm 



'^2 01 



< /(r) < 1 
7 



for all r. 



(23) 



Our focus in this paper will be hypothetical continuous functions /(r) that depend on the 
distance r = |r| and could potentially correspond to statistically homogeneous and isotropic 
media without long-range order such that /(O) = 1 and /(r) tends to zero as r — oo sufficiently 
fast so that the Fourier transform of x{f^) = *S'2*^(r) — (pf exists. The latter two properties of 
/(r) ensure that 5*2*^ (r) obeys its proper asymptotic limiting behaviors as specified by (jH)) and 
(ffj). respectively. When the scaled autocovariance /(r) depends only on the magnitude r = |r|, 
then the Fourier transform condition (jSJ on f{k) can be written in any space dimension d 

as [221 

m = {2.)i [r'^-^f{ry^^^dr>^^ (24) 
where /c = |k| and Jy{x) is the Bessel function of order v. The bounds are equivalent to 



mm 



')2 01 



< Unf < 0, (25) 



where is the infimum of /(r). Note that when function /(r) is independent of the volume 
fraction 0i, it would correspond, if realizable, to a two-phase medium with phase-inversion 
symmetry. 3j A two-phase random medium possesses phase-inversion symmetry if the geome- 
try of phase 1 at volume fraction 0i is statistically identical to that of phase 2 in the system 
where the volume fraction of phase 1 is 02 and hence 

5(^)(r,0i,02) = 5f (r;02,0i). (26) 

By construction, the upper bound of is always satisfied. All of the functions /(r) consid- 
ered below are taken to be independent of the volume fraction 0i and therefore any violation 
of the lower bound of (j^^ implies that a two-phase statistically homogeneous and isotropic 
medium cannot exist for the following volume-fraction intervals 

O<0.< , i^i"/' I and ^ ^ I „ I < 0. < 1. (27) 

Figure 121 depicts the bounds on for /(r) that are independent of volume fraction. 

First we note that for any /(r) that monotonically decreases in r to its long-range value 
of zero, the pointwise nonnegativity condition is obeyed for < 0i < 1. However, as 
some examples below will demonstrate, such an /(r) does not necessarily obey the triangular 
inequality ([221) • ^ natural example of a monotonic scaled autocovariance function /(r) is the 
simple exponentially decaying function, i.e., 

/(r) = exp(-r/a), (28) 
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Figure 2: The bounds ()25|) on the infimum for volume-fraction- independent f{r). The 
quantity finf must he within the region dehneated by the heavy sohd hues. 

where a is a positive parameter that we caU the "correlation length." This function was first 
proposed by Debye and coworkers, [TTJ who intuited that it should correspond to structures 
in which one phase consists of "random shapes and sizes," but presented no proof that such 
was the case. The function ()28|) obeys the necessary nonnegativity condition on the 
spectral function f{k) for any d as well as the triangular inequality (f^ . The satisfaction of 
these necessary conditions does not ensure that such a correlation is realizable. However, the 
aforementioned inverse optimization construction technique [7| |3] was applied to generate a 
two-dimensional digitized realization corresponding to (PHj) (see Figure E)). This leads one to 
believe that ()28|) is exactly realizable. Indeed, there are specific two-phase microstructures 
that achieve the "Debye" random-medium function (PHj) in the plane. The function (j^H|) 
is a special case of a more general realizable subclass of B given by the completely monotonic 
functions,[Wi i-e.. 



where F{X) is a nonnegative bounded measure [bounded and nonincreasing function on 
(0, oo)], i.e., dF>0 and /o°° dF{\) = 1. We see that if F = e(A - a"^), then dF = S{\ - a'^) 
and ()28|) is recovered, where G(x) and 6{x) are the Heaviside and Dirac delta functions, 
respectively. 

Another natural monotonic scaled autocovariance function /(r) to consider is the Gaussian 
function, i.e.. 




(29) 



/(r) = exp[-(r/a)^]. 



(30) 



Although any such Gaussian function has a nonnegative spectral function f{k), it cannot 
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Figure 3: Construction of a digitized two-dimensional realization of a "Debye" random medium 
(400 X 400 pixels) . [3 Ej Here the volume fraction 0i = 02 = 0.5 and correlation length a = 2 
pixels. 



correspond to a two-phase random medium in because the slope of (r) at r = is zero 
(i.e., the specific surface is zero) and therefore violates condition (fT^ or, more generally, the 
triangular inequality (j22I)- For precisely the same reasons, the class of monotonic functions 



/(r) = exp[— (r/a)"]. for any a > 1 (31) 

and 

^^""^ " [1 + (r/a)^]/^-i any /3 > d (32) 

cannot correspond to a two-phase random medium in d dimensions. These specific examples, 
some of which are illustrated in Fig. |5 show that the nonnegativity condition ()24|) and 
triangular inequality (fT^ are independent necessary conditions. 

The final monotone function that we test is the simple linear function 

{1 — r/a, if r < a, 
' ' - ' (33) 

0, otherwise, 

Shepp (IHl proved that such a scaled autocovariance is realizable by a statistically homogeneous 
two-phase medium in one dimension. However, this autocovariance is not realizable in higher 
dimensions because its spectral function f{k) can take on negative values for certain values of 
k. It is noteworthy that it has been shown that for any positive definite /(r) in one dimension, 
the function 2 arcsin(/)/7r as well as Stt^^ EkC^k + 1)V((2A; + l)r) are in B. 

A generalization of the Debye random-medium function ()28|) that is nonmonotone and 
would be characterized by short-range order is the following expression: ||9] 

;(,) = e-'" (34) 
qr 
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Figure 4: Examples of scaled autocovariance functions that cannot correspond to statistically 
homogeneous and isotropic two-phase random media: (1) /(r) = 1/(1 + r"^)^ for d < 5; (2) 
/(r) = exp(— r^) for any d; and (3) /(r) = exp(— r®) for any d. 

where q is an inverse length scale that controls oscillations in the term sin(gr)/(gr). The 
spectral function f{k) of (j34|) in one, two, and three dimensions obeys the nonnegativity 
condition (j^ . Interestingly, Torquato [3] observed that although (jMj) satisfies the upper 
bound of binary condition (fTTj) . it does not necessarily satisfy the lower bound of (fTTj) or, 
equivalently, the lower bound of for all 0i, depending on the values of a and q. In other 
words, there are values of the infimum which in this case is a true global minimum, that 
violate the lower bound of ()25|) . Let tq be the radial distance at which /(r) achieves its global 
minimum. The minima of /(r) are solutions to the transcendental equation q{a + r) tan(gr) = 
g^ar. The extremum value qr^ can be shown to lie in the interval [7r,37r/2) for arbitrary a 
and q. For example, for aq = Svr, vq ^ 5.671a and fmin ~ —0.1818 (see Fig. For example, 
for aq = 8tt, vq ^ 0.1772a and fmin ~ —0.1818 (see Fig. El), and therefore, according to (j2Z|), 
(jH^ is not realizable for the volume-fraction intervals 

< 0i < 0.1538... and 0.8461 ...< 0^ < 1. (35) 

Interestingly, two realizations of digitized two-dimensional two-phase media were previously 
constructed EI that putatively correspond to the scaled autocovariance function (jMj) for 
02 = 0.2 and 0.5, respectively, and the aforementioned choice of a and q are shown in Figure IHl 
At 02 = 0.2, the system resembles a dilute particle suspension with "particle" diameters 
of order b. At 02 = 0.5, the resulting pattern is labyrinthine such that the characteristic 
sizes of the "patches" and "walls" are of order a and 2n/q, respectively. For these sets of 
parameters, all of the aforementioned necessary conditions on the function are met, except 
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0.2 0.4 0.6 0.8 1 

r/a 

Figure 5: The damped sinusoidal function ()34|) with qa = Svr. 

for the triangular inequality. Although (jMjl satisfies the negative slope condition (fT^ at the 
origin, it only satisfies the convexity condition ()17|) for qa < -\/3, which we see is violated in 
these instances, implying that the triangular inequality must be violated. As it turned out, 
the construction procedure matched the target function (jMj) for almost all r, but it could not 
yield convex behavior in the vicinity of the origin. Since the triangular inequality was not 
known at the time, it was difficult to ascertain whether the slight discrepancy in the curvature 
of the function at the origin was numerical imprecision. We now know in retrospect that the 
construction technique was revealing that a two-phase medium with a scaled autocovariance 
function (jH^ cannot be exactly realized, which is a testament to the power of this method. 

4 Conclusions 

We have identified all of the known necessary conditions on the two-point correlation func- 
tion 5*2*'' (r) of statistically homogeneous two-phase media and have derived a new but simple 
positivity condition that it must satisfy. Using these conditions, we were able to ascertain 
the realizability of certain classes of proposed correlation functions. In future work, it will be 
important to identify other checkable necessary conditions. The stochastic optimization con- 
struction technique appears to be a very powerful numerical tool in guiding such a search. 
Finally, we note that the analogous realizability problem for the pair correlation function g2 of 
point processes |13| I14[ IT^ I24[ 1^ offers many interesting challenges. It has recently been con- 
jectured that the known standard non-negativity conditions on g2 are sufficient to ensure the 
existence of point processes at and above some sufficiently high space dimension [211 12ZI- Appli- 
cation of this conjecture implies the possibility that the densest sphere packings in sufficiently 
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02 = 0.2 (t)2 = 0.5 

Figure 6: Construction of digitized two-dimensional realizations (400 x 400 pixels) that 
putatively correspond to the target function given by ()H4|1 for 02 = 0.2 and 0.5.|nilSl Here 
a = 32 pixels and g = Svr/a. We now know that this function is not exactly realizable because 
even though the construction technique matched ()34p for almost all r, it could not yield the 
necessary convex behavior in the vicinity of the origin. 

high dimensions are disordered rather than periodic, implying the existence of disordered clas- 
sical ground states for some continuous potentials. In future work, it would be interesting to 
investigate whether an analogous conjecture applies to binary stochastic processes. 
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